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We discuss a piecewise-conserved constant of motion for a simple dissipative oscillatory mechani- 
cal system. The system is a harmonic oscillator with sliding (dry) friction. The piecewise-conserved 
constant of motion corresponds to the time average of the energy of the system during one half- 
period of the motion, and changes abruptly at the turning points of the motion. At late times the 
piecewise-conserved constant of motion degenerates to a constant of motion in the usual sense. 

I. INTRODUCTION 

Finding constants of motion is an important step in the solution of many problems of Physics, as they allow to reduce 
the number of the degrees of freedom of the problem. Constants of motion are intimately related to conservation 
laws or symmetries of the system. For example, it is well known that a symmetry of the Lagrangian of a system 
is responsible, by virtue of Noether's theorem, to a constant of motion [Q. By definition, a constant of motion 
preserves its value during the evolution of the system. Even in cases where there are no constants of motion, one 
could still sometimes find adiabatic invariants, which generalise the concept of a constant of motion to systems with 
slowly-varying parameters. It turns out, however, that there are cases where there are constants of motion, which 
are only piecewise-conserved. Well-known examples of such piecewise-conserved constants of motion arise from the 
generalisation of the Laplace-Runge-Lenz vector for general central potentials. For example, in the case of the three- 
dimensional isotropic harmonic oscillator, the Fradkin vector ^ (which generalises the Laplace-Runge-Lenz vector) 
abruptly reverses its direction (although preserves its magnitude) during a full period H (see also Q|). (The Fradkin 
vector directs toward the perigee, and the position of the perigee jumps discontinuously whenever the particle passes 
through the apogee.) Also, in the truncated Kepler problem, the Peres- Serebrennikov-Shabad vector abruptly changes 
its direction whenever the particle in motion passes through the periastron Again, it is just the direction of the 

conserved vector which is only piecewise-conserved: the magnitude of the vector remains a constant of motion in the 
original meaning (namely, the magnitude has a fixed value throughout the motion). Piecewise-conserved constants 
of motion may also be relevant for systems which involve radiation reaction. In the above examples, the piecewise- 
conserved constants of motion appear in non-dissipative systems, and result from a discontinuity of the force (as in 
the truncated Kepler problem) or from geometrical considerations (as in the three-dimensional isotropic harmonic 
oscillator case). These piecewise-conserved constants of motion involve vectors rather than scalars, which still conserve 
their magnitude. 

The following question arises: Can one find, for elementary systems, piecewise-conserved constants of motion? In 
what follows, we shall discuss an elementary piecewise-conserved scalar constant of motion, for a simple oscillatory 
mechanical model which involves dissipation in the form of sliding (dry) friction. Although dry friction is more nearly 
descriptive of everyday macroscopic motion in inviscid media, it is usually ignored in elementary mechanics courses 
and textbooks, which very frequently discuss viscous friction (which is velocity dependent). Dry friction exhibits, 
however, some very interesting features and can be readily presented in the laboratory or classroom demonstration. 

The problem of a harmonic oscillator damped by dry friction was considered by several authors. Lapidus analised 
this problem for equal coefficients of static and kinetic friction, and found the position where the oscillator comes to 
rest [Q. Hudson and Finfgeld were able to find the general solution of the equation of motion M. However, they again 
assumed equal coefficients of static and kinetic friction, and used the Laplace transform technique, which is unknown 
to students of elementary mechanics courses, to generate the solution. An elementary solution, which ignores static 
friction, was derived by Barratt and Strobel || . This solution is based on solving separately for each half cycle of the 
motion, and is consequently tedious and unappealing. Recently, Zonetti et al. [[l0| considered the related problem of 
both dry and viscous friction for a pendulum, but did not offer a full analytic solution for the motion. 

In this work, we find the general solution for the motion taking into account both static and kinetic friction, using 
elementary techniques which are available for students of elementary courses of mechanics. We analyse the solution 
using a piece-wise conserved constant of motion. The discussion, as well as the corresponding laboratory experiment 
of classroom demonstration, are suitable for a basic course for physics or engineering students. 
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II. ELEMENTARY DISCUSSION 



Let a block of mass M be placed on a horizontal surface, such that the coefficients of static and kinetic friction are 
ji s and /ik, respectively. The block is attached to a linear spring with spring constant k (for both compression and 
extension), such that initially the spring is stretched from its equilibrium length by £, and the block is kept at rest 
at x = I. At time t = the block is released. If k£ > fi s Mg, g being the gravitational acceleration, the block would 
start to accelerate. We take the friction to be small, namely fi s ,k -C £k/(Mg), and also assume slow motion, such 
that the friction force is independent of the speed. Namely, we neglect effects such as air-resistance, and include only 
the force which results from the block touching the surface. We also neglect any variation of Hk with the speed. 

Immediately after the block starts accelerating, its motion is governed by the equation of motion 

Mx = -kx + fi k Mg, (1) 

with initial conditions a;(0) = £ and i(0) = 0. From now on, let us introduce the frequency lo 2 — k/M . Of course, the 
system does not preserve its energy, due to the friction force. However, let us define a new coordinate x' — x — [i^g j io 1 . 
Equation (p then becomes x' + uj 2 x' = 0. For this equation we know that there is a constant of motion, namely 

£ = ^Mx' + ^Mcu 2 x' 2 . Therefore, despite the presence of friction, one can still find a constant of motion, which 
has the functional form of the total mechanical energy, but which is of course not the energy, as the latter is not 
conserved. Calculating its numerical value we find that £q = |Mw 2 {£ — pu^gjuj 2 ^ . 

At the time t = ir/w the velocity of the block vanishes, and it can be easily shown that at t — ir/w its acceleration 
is x(tt/uj) = (£ — fikg/ijJ 2 )^ 2 > 0, such that the block reverses its motion. (We assume here that Mx(tt/uj) > fi s Mg.) 
The nature of the friction force is that its direction is always opposite to the direction of motion. Consequently, the 
equation of motion now changes to 

Mi = -kx - n k Al g, (2) 

with initial conditions x(ir/tu) = £ — 2[ikg/u 2 and x(tt/uj) = 0. One can again solve this equation readily. This time, 
let us define x' = x + iikg/io 2 . Equation (||) again becomes x' + ui 2 x' = 0, such that £ is still conserved. However, 
this time the numerical value of £, which we denote by £\ ^ £q, and we find £\ — ^Mu> 2 {£ — 3fikg/^J 2 ) ■ One can 
describe the next phases of the motion similarly. During each phase of the motion (during half a period between two 
times at which the velocity vanishes) £, if defined properly, is conserved. However, £ is only piecewise-conserved, 
as its value changes abruptly from phase to phase. We note that the period of the oscillations is not altered by the 
presence of friction, and denote by P\/2 half that period. Namely, P 1 / 2 = it/uj. 



III. GENERAL DISCUSSION 



Let us now discuss the system in a more general way. It turns out that although there are friction forces, one can 
still write a hamiltonian 

H{x,p-t) = ^p 2 + X -Muj 2 x 2 - f(t)x, 

where f(t) = (— l)^ Pl/2 l(iiMj. The equation of motion is now 

x + uj 2 x = f(t)/M (3) 

with the initial conditions being (as before) a;(0) = £ and i(0) = 0. We denote by square brackets of some argument 
the largest integer smaller than or equal to the argument. We also assume that the static friction force at the turning 
points of the motion is smaller than the elastic force of the spring, such that the motion does not stop. (Of course, for 
large enough time, this would not be true any more, and the block would eventually stop — see below.) Let us define 
the (complex) variable £ = x + icux where i 2 = —1. Then, instead of a real second order equation (such as Eq. 

), one obtains a complex first order equation. It is advantageous to do this, because there is a general solution for 
any inhomogeneous linear first-order differential equation in terms of quadratures. Substituting the definition for £ 
in Eq. (||) we find that the equation of motion, in terms of £, takes the form 

e - *we = m/M, (4) 
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with the initial condition £o = £(t = 0) = iiol. The solution of the equation of motion (0) is 

m = ^ {Jo dt ' + C °} ' (5) 

After finding the solution we can find x(t) and x(t) by x{t) = (£(*) + £*(*))/ 2 and x(t) = (£(t) - £*(t))/(2i). We 
denote by a star complex conjugation. In order to integrate Eq. (|3|) we find it convenient to separate the discussion 
to two cases: case (a) where [t/Pi/2] is an odd number (namely [i/Pi/g] = 2n — 1), and case (b) where [V-P1/2] 
is even (namely, [t/Pi/2] — 2n), where n is integer. We next split the interval of integration in Eq. (^J) into two 
parts: we first integrate from t' = until t 2n -i = (2n — l)P 1 / 2 , and then integrate from t2n-i to t, and sum the two 
contributions. Integrating term by term we find 

'.2n-i 1 r (2j+l)P 1/2 



-/(i')e- Mt dt' = J2 Hk9e- iut dt' 

M ~" J2jP 1/2 



n-2 „(2j+2)P 1/2 



„-_n J(2l + l)Pi 



For case (a) we find that 



J=0 ^(2i+l)Pi /2 

—2i(2n—l)(i k g/<-^- (6) 



■!/(;>— *' dt' = - / /i fc ge-^ 4 ' dt' - -i(e- iut + l)n k g/w. (7) 



t a »- 1 M Jt^-i 



For case (b) we find 

-J(t')e- iut 'dt' = -l "' fJ, k ge~ iut ' dt' + I n k ge~ iwt ' dt' 



-1 p2nP 1 / 2 rt 



a«,-l J(2n-l)P 1/2 J2nP 1/2 

= i(e- 1 ^ - 3)^ k g/u. (8) 
Collecting the two integrals, we find for case (a) that 

£„(t) = -2i(2n - l^Ukg/u - i(l + e iwt )n k g/oj + iwfo 4 "*, (9) 

and for case (b) 

&(t) = -2i(2ra - l)e^Vfe.9/^ + *(1 - 3e ia, *)/x & 3/w + ^e iwt . (10) 

Recalling the different values of [t/Pi/ 2 ] for the two cases (a) and (b), we can unify the expressions for both £ a (f) 
and £,b(t), namely 

£(t) = -i (2 [t/P 1/2 ] + 1) e iut ii k g/w + (-l)i t/Pl ^i^ k g/uj + iujle wt . (11) 
From this solution for £(i) we can find that 

x(t) = (-lp p ^}fi k g/u) 2 + {£ - (2 [t/Pi/a] + l) /W^ 2 } cosc^t, (12) 

and 

i(t) = (2 [t/Pi/ 2 ] + l) Mfc.g/w 2 } wsinwt. (13) 



An interesting property of the solution given by Eqs. (|12|) and ( |13[ ) is that for each half cycle it looks as if the motion 
were that of a simple harmonic oscillator, with no friction. In fact, the effect of the friction for each half cycle enters 
only in the initial conditions for that half cycle, or, more accurately, in the smaller value for the initial position for 
the half cycle. In addition, it is evident from Eq. ( |l2| ) that the damping of the amplitude of the oscillation is linear 
in the time t, whereas in disspative systems in which the resistance is speed-dependent the damping is exponential in 
the time. 
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Let us now define a new coordinate x'(t) = x(t) — f(t)/(Mui 2 ) = x(t) — (— l^^^ifi^g/uj 2 . Then, we find that 

x'{t) = {£ - (2 [t/P 1/2 ] + 1) Hkg/u 2 } coscut, (14) 

and 

x>(t) = - {£ - (2 [t/Pi/a] + l)^ k g/uj 2 }ujsmujt. (15) 



Next, we define 



£{t) = ^Mx ,2 (t) + ^Muj 2 x' 2 (t). (16) 



Substituting the expressions for x'(t) and x'(t) in £, we find that 

£{t) = X -Muj 2 {£ - (2 [t/P 1/2 ] + 1) fiug/u: 2 } 2 . (17) 

It is clear that £ is not a constant of motion. However, a close examination shows that it is piecewise conserved: the 
only dependence on t is through [i/P^]- Therefore, between any two consecutive turning points we find that the 
numerical value of £ is conserved. Consequently, £ is a piecewise-conserved constant of motion. 

Clearly, £ has the dimensions of energy. However, we stress that £ is not the mechanical energy of system, 
because the latter is not even piecewise-conserved. In fact, the total mechanical energy of the system is T{t) — 
\Mx 2 (t) + \Muj 2 x 2 {t), namely, 

T(t) = 1 -Muj 2 {£ (2 [</P 1/2 ] + 1) » k g/to 2 } 2 + \u^f 

A A UJ 

+ (-ift/FiM^Mg {£ - (2 [t/Pi/a] + 1) /W^ 2 } coso;i, (18) 

which is a monotonically decreasing function of t, as expected. (Notice that whenever the cosine changes its sign, so 
does its amplitude.) Of course, if we add to T(t) the work done by the friction force, we obtain a constant value. The 
fact that the total mechanical energy is monotonically decreasing is important: the system loses energy constantly. 
We have previously noted that the position and the velocity of the block during each half cycle are influenced by 
the presence of friction only through the initial conditions for that half cycle, but otherwise the motion is simple 
oscillatory. Despite this fact, the loss of energy occurs throughout of motion, as is evident from Eq. (|l8|), as should 
be expected. 

In order to gain some more insight into the meaning of the piecewise-conserved £ , let us find the time average of 
T(t) between two successive turning points. Clearly, the average of the cosine vanishes, and we find 

< T{t) > = \Mu? {£ - (2 [t/P 1/2 ] + 1) ^g/u 2 } 2 + \m^- 

A A LU 



£{t) + \M^f. (19) 



Therefore, the physical meaning of £ is the following: up to a global additive constant (namely, a constant throughout 
the motion) £ is equal to the time average of the total mechanical energy of the system T(t) between any two 
consecutive turning points. Because of the dissipation, this time average decreases from one phase of the motion to 
the next, and therefore £ is only piecewise conserved. 

We next present our results graphically for two sets of parameters. First, we choose the parameters £ — 1 m, 
uj = 5scc _1 , M = 1 Kg, g = 9.8m/sec 2 , fi s = 0.54, and fi^ — 0.36. (These values for the coefficients of friction 
are typical for copper on steel.) In all the figures below the units of all axes are SI units. Figure [l] displays the 
position x(t) and the velocity x(t) vs. the time t. It is clear that the amplitude of the oscillation attenuates, and 
eventually the block stops in a state of rest. Figure § displays the piecewise-conserved £ and the mechanical energy 
T as functions of the time t. Indeed, the energy T(t) is a monotonically-decreasing function of t, whereas £(t) is 
piecewise-conserved. One can also observe that up to a constant indeed £ is the average of the energy T over one 
half-cycle of the motion. The dissipation of energy is most clearly portrayed by means of the phase space. Figure || 
shows the orbit of the system in phase space, namely, the momentum p — Mx vs. the position x. The loss of energy 
is evident from the inspiral of the orbit. Eventually, the orbit arrives at a final position in phase space, and stays 
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there forever. For figures UJ|, and ^ we changed only the coefficients of friction to (i s = 0.15 and fXk = 0.06. (These 
parameters are typical for the contact of two lubricated metal surfaces.) As the coefficients of friction in this case 
are smaller than their counterparts in the former case, we can observe many more cycles of motion before the motion 
stops. (In fact, the number of half-cycles in this case agrees with Eq. (|o|) below.) We note that because of the scale 
of Fig. U it is not apparent that the energy arrives at a non-zero constant value at late times. In this case, also, 
the qualitative characteristics of the motion are the same as in the former case (Figs. but here the attenuated 

oscillatory motion is more apparent. 




t 

FIG. 1. The position and the velocity as functions of time t, for the choice of parameters fj, s 
line: x(t). Dashed line: x(t). 



0.54 and fi h = 0.36. Solid 



f 

FIG. 2. The picewise-conserved constant of motion and the total mechanical energy as functions of time t, for the choice of 
parameters fj, a = 0.54 and /i^ = 0.36. Solid line: £ (t). Dashed line: T(t). 
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FIG. 3. The orbit of the system in phase space: the momentum p = Mx vs. the position x, for the choice of parameters 
fi s = 0.54 and (j,k = 0.36. 
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FIG. 4. The position and the velocity as functions of time t, for the choice of parameters fi s — 0.15 and fi k — 0.06. Solid 
line: x(t). Dashed line: v = x(t). 
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FIG. 5. The picewise-conserved constant of motion and the total mechanical energy as functions of time t, for the choice of 
parameters fj, s = 0.15 and fj,k = 0.06. Solid line: £ (t). Dashed line: T(t). 
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FIG. 6. The orbit of the system in phase space: the momentum p = Mx vs. the position x, for the choice of parameters 
fj, s = 0.15 and jik = 0.06. 



Of course, the motion will not continue forever: Because of the decrease in the amplitude of the motion, eventually 
the static friction force at some turning point would be larger than the elastic force exerted on the block by the spring. 
Namely, at t = nP X j 2 we find x(t = nPi/ 2 ) = (— l) n (£ — 2nfj,kg /u) 2 ) , for some integer n, and the motion will stop for 
fi s g > ui 2 (£ — Innkd /w 2 ) , or after an integral number of phases which is equal to the least integer n which satisfies 



-> 1 ( uH 
n > - 

2 VMfcff 



Ms 



(20) 



We note that for the special case where oj 2 £ 2 / (nug) is integral the block may stop at x = 0. This happens, however, 
only for special values of the parameters of the systems, and in general the system will rest at x ^ 0. 

Then, the block would remain at rest, and £ would be a constant of motion from then on. Namely, because of the 
dissipative nature of the problem, eventually the piecewise-conserved constant of motion becomes a true constant of 
motion, but this happens only when the dynamics of the system becomes trivial. (In our case, when the system is in a 
constant state of rest.) This feature of the dissipative system is in contrast with other piecewise-conserved constants 
of motion, which arise from non-dissipative systems, such as the truncated Kepler problem or the three-dimensional 
isotropic harmonic oscillator, where the piecewise-conserved constant vector remains piecewise conserved for all times. 
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